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Abstract. We show how the interplay of non-linear dy- 
namics, self-gravity, and fluctuations leads to self-affinc 
behavior of matter density correlations quite generically, 
i.e., with a power law exponent whose value does not de- 
pend in a very direct way on the specific details of the 
dynamics. This we do by means of the Renormalization 
Group, a powerful analytical tool for extracting asymp- 
totic behavior of many-body systems. 

Key words: Gravitation - Hydrodynamics - Methods: 
analytical - Cosmology: theory - Large-scale structure of 
Universe 



to condensed matter physics (Wilson & Kogut 1974), later 
on extended to problems in hydrodynamics (Forster, Nel- 
son and Stephen 1977) and applied to surface growth 



1. Introduction. 

The Universe at short and moderate distance scales is in- 
homogeneous, being filled by numerous structures. In fact 
the clumpy structure of matter in the Universe extends for 
over 15 orders of magnitude in linear size, from stars to 
the largest clusters of galaxies, and for about 18 orders of 
magnitude in mass (Ostriker 1991). At scales from galax- 
ies to clusters of galaxies and even superclusters, the Uni- 
verse exhibits self-similar, fractal behavior (Mandelbrot 
1983; Peebles 1993). At the largest scales, those probed 
by COBE observations (Smoot 1991), for example, the 
Friedmann-Robertson- Walker (FRW) cosmology which is 
based on the hypotheses of homogeneity and isotropy, pro- 
vides adequate description of many uncorrelated observa- 
tions. Isotropy is broken by the primordial fluctuations in 
the matter density which subsequently become amplified 
by the gravitational instabilities. 

In this paper we explore the observed self-affinity and 
scaling in the density correlation function and tackle them 
with the renormalization group, a tool originating in quan- 
tum field theory (Gell-Mann & Low 1954), then extended 



(Kardar, Parisi and Zhang 1986) and in the last few years 
ap plied to gra vitation and cosmology (Perez-Mercader et 
al. |1996, 19971 ) . 
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Our starting point will be the non-relativistic hydro- 
dynamic equations governing the dynamical evolution of 
an ideal self-gravitating Newtonian fluid in a FRW back- 
ground. Through a series of steps involving physically jus- 
tified approximations, these equations reduce to a single 
evolution equation for the cosmic fluid's velocity potential. 
At this point one has a single deterministic hydrodynamic 
equation. To this we add a noise source which represents 
the influence of fluctuations and dissipative processes on 
the evolution of the fluid, arising from, but not limited 
to, viscosity, turbulence, explosions, late Universe phase 
transitions, gravitational waves, etc. The resulting equa- 
tion is recognized as a cosmological variant of the K ardar- 
Parisi-Zhang (KPZ) equation (Kardar et al. 1986 ) which 
is the simplest archetype of nonlinear structure evolution 
and has been studied extensively in recent years in the 
co ntext of surface growth phenomena (Barabasi & Stan- 
ley 1995 ). While the essential steps involved in arriving at 
this KPZ equation are reviewed in Section ^, we refer the 
reader to the concise details of its complete cosm ological- 
hydrodynamical derivation in Buchert et al. |1997 . The rel- 
evance of stochastic fluctuations in structure formation in 
the Universe has been addressed in Berera & Fang 1994 . 
where the role of the KPZ equation was emphasized. For 
FRW cosmologies with curved spatial sections, this equa- 
tion contains a time-dependent mass-term. In Section || we 
discuss briefly the adiabatic approximation which allows 
one to treat this mass term as a constant during most of 
the expansion of the background cosmology. 

Having thus reduced the essential hydrodynamics to a 
single dynamical stochastic equation, we turn to the re- 
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suits of a dynamical renormalization group analysis to cal- 
culate the precise form of the density-density correlation 
function using the KPZ equation. To do so we proceed 
first by extracting the renormalization group equations 
governing the evolution of couplings appearing in the cos- 
mological KPZ equation with respect to changes in scale 
and integrating out the short distance degrees of freedom. 
As we are interested in the behavior of the system at ever 
larger scales, we will run the renormalization group (RG) 
equations into the large scale (infrared) limit, a process 
known as coarse-graining. The fixed points of these RG 
equations determine the long-time, long-distance behav- 
ior of the system and their subsequent analysis yields both 
the power law form of the correlation functions as well as 
their explicit exponents. In particular, the value of the 
galaxy-galaxy correlation function exponent 7 ~ 1.8 is 
calculated, and its value understood in terms of spatio- 
temporal correlated noise. A summary of our results and 
discussion are presented in the closing section. 

2. Prom Hydrodynamics to the KPZ Equation. 

Structure formation in the Universe in the range of a few 
Mpc up to several hundred Mpc can be modelled as the 
dynamical evolution of an ideal self-gravitating fluid and 
studied within the framework of General Relativity. How- 
ever, we will restrict our treatment to the case of New- 
tonian gravity for a number of physically justifiable rea- 
sons. First and foremost, the length scales involved in large 
scale structure formation by matter after decoupling from 
radiation are smaller than the Hubble radius (of order 
w 3000Mpc at the present epoch), so that general rel- 
ativistic effects are negligible. Secondly, it is known ob- 
servationally that matter flow velocities are much smaller 
than the speed of light and that this non-relativistic mat- 
ter plays a dominant role over radiation in the epoch of 
structure formation. Thus, the evolution in space and time 
of large scale structure in the Universe will be described 
adequately by the Newtonian hydrodynamic equations for 
a fluid whose components are in gravitational interaction 
in an expanding background metric []. The connection with 
observations and the subsequent phenomenology is made 
by writing these equations in terms of the comoving coor- 
dinates x and coordinate time t, and are (Peebles 1992) 



V • w = -4irGap(t)5, V x w = 0, 



(3) 
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1 The resolution length must be large enough so to as to 
legitimize the use of continuum equations (i.e., to validate the 
long-wavelength, or hydrodynamic limit) for the description of 
the matter density, peculiar velocity and acceleration fields. 
This requirement implicitly singles out a lower length scale 
which we may take to be greater than roughly the mean galaxy- 
galaxy separation length (~ 5 Mpc). 



where v(x,t) and w(x,t) are the peculiar velocity and 
acceleration fields of the fluid, respectively. In the above 
equations, H = a(t)/a(t) is the Hubble parameter, a(t) 
is the scale factor for the underlying cosmological back- 
ground and p(t) is the (time-dependent) average cosmo- 
logical density. Thus, the local density is p — p(l + S) and 
the dimensionless density contrast is 8 = 

In the linear regime and for dust (p — 0), the vortic- 
ity of the velocity field v damps out rapidly by virtue of 
the background expansion, and therefore v can be derived 
from a velocity potential ip, v — — Vf/', in the long-time 
limit. Furthermore, parallelism between w and v holds, i. 



w = F(t)v 



(4) 



where F(t) is the long-time limit solution of the following 
Ricatti equation 



F = AnGp{t) - HF - F 2 



(5) 



Following the Zel'dovich approximation, we make the 
assumption that in the weakly non-linear regime, the par- 
allelism condition (|]) continues to hold, with F(t) given as 
the solution of (pb but we now employ the fully non-linear 
Euler equation (0) with pressure term included. That is, 
we assume that the non-linearities and pressure have not 
yet had enough time to destroy the alignment in the accel- 
eration and velocity fields of the cosmic fluid. Assuming 
also that the pressure is a function of the density, p = p{p) , 
one arrives at the following equation for the velocity po- 
tential (Buchert et al. |1997|) 



^ = vfi(t)V 2 ip + ~/ 2 (i)(W) 2 + + ft(x, t). (6) 



A 
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T 



where we have also added a stochastic source or noise r] 
(see a few lines below for a discussion). The three dimen- 
sionless functions of time, fx, / 2 and / 3 are 
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p'(p)F(t) 1 
AnGp(t)a 2 (t) v ' 

1 



a(t)y 
f 3 (t) = (F(t) - H{t))T. 



(7) 

(8) 
(9) 



The positive constants A and T are introduced to 
carry dimensions. One can think of them as typical val- 
ues of the corresponding time-dependent coefficients dur- 
ing the epoch we are interested in. The reason why we 
introduce them here explicitly (instead of setting them 
equal to unity) will become clear in Section [|. The f\- 
term arises from the pressure term in the Euler equation 
which we have expanded to lowest order about the zero- 
pressure limit. Note that higher-order terms in this Taylor 
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expansion will yield higher-derivative terms (0(V 2 ip) 2 ) in- 
volving quadratic and higher powers of the field ip. The 
/2-term is simply the convective term from the original 
Euler equation (Q) written in the new variables, while the 
/3-term entails the competition between the damping of 
perturbations due to the expansion (—Hip) and the en- 
hancement due to self-gravity (+Fi/j). 

The Zel'dovich approximation iQ) together with the 
Poisson equation yield an important relationship between 
the density contrast and V 2 ip, the Laplacian of the veloc- 
ity potential, namely 



d(x,i) 



F 



AirGap 



(10) 



which shows that the density contrast tracks the diver- 
gence of the peculiar velocity. Because of this identity, we 
are able to calculate density correlation functions in terms 
of (derivatives of) velocity potential correlations. This is 
why it is worthwhile investigating the scaling properties 
of the KPZ equation in detail. 

The noise rj(x,t) appearing in (||) represents the effects 
of random forces acting on the fluid particles, including the 
presence of dynamical friction, and also of degrees of free- 
dom whose size is smaller than the coarse-graining length 
(indeed, the very fact that we are using continuous fields 
to describe the dynamics of a system of discrete particles 
means that we are implicitly introducing a coarse-graining 
length, such that the details below this resolution length 
are not resolvable. (The relevance of this length-scale will 
become clear in section ^j.) There exist a number of phys- 
ical processes on various length and time scales that con- 
tribute to an effective stochastic source in the Euler equa- 
tion. Indeed, any dissipative or frictional process leads to 
a stochastic force, by virtue of the fluctuation-dissipation 
theorem. So, fluid viscosity and turbulence should be ac- 
countable to some degree by adding a noise term in the 
dynamical equations. Early and late Universe phase tran- 
sitions, the formation of cosmic defects such as strings, 
domain walls, textures, are sources for a noisy fluctuating 
background, as are also the primordial gravitational waves 
and gravitational waves produced during supernovae ex- 
plosions and collapse of binary systems. Another way to 
visualize the noise is to imagine a flow of water in a river; 
there are stones and boulders in the river-bed which obvi- 
ously perturb the flow and have an impact on its nature. 
The noise can be thought of as a means for modelling the 
distribution of obstacles in the river bed. In cosmology the 
obstacles in the river-bed represent, e.g., density fluctua- 
tions and the fluid is the matter flowing in this "river" . 
The noise is phenomenologically characterized by its av- 
erage value and two-point correlation function as follows, 



( V (x,t))=0 

(r)(x,t)r)(x',t')) = 2D(x,x';t,t'), 



(11) 
(12) 



with D{x,x';t,t') a given function of its arguments (see 
below). All higher cumulants vanish. The noise is thus 
Gaussian. However, the velocity field need not be (and in 
general, will not be) Gaussian as a consequence of coarse- 
graining of the dynamics. 

The stochastic hydrodynamic equation @ can be fur- 
ther cast into 

^ = zA7 2 * + ^A(W) 2 - to 2 (t)* + rj(x, r), (13) 

by means of the following simple change of variables and 
redefinition of the "physical-time" t into a "conformal- 
time" t 



^( x ,r)^^f^Mx,t(r)), 



= / dt'h(t'). 



The linear (mass) term in equation (|T^) is given by 



m 2 (r) = In 

or 



Mr) / 1 

ex P 



h(r) 



T 



t(r) 



dt' h{t') 



to (to) 



and the noise transforms into 

fj(x,t) = M±r,(x,t). 
hit) 



(14) 
(15) 

(16) 
(17) 



Finally, Eq. (|lj) is a generalization to a cosmological 
setting of the Kardar-Parisi-Zhang (KPZ) equation for 
surface growth (Kardar et al. 1986|) , and differs from the 
standard KPZ equation in the linear (mass) term which, 
as seen from Eq. jl6|), originates in the expansion present 
in the background cosmology. Here, v plays the role of a 
diffusion constant while A is proportional to the average 
"speed of growth" of ip, 

3. Adiabatic Approximation. 

In this section we briefly comment on some results about 
the time-dependent mass term (|l^) of equation ([[3]) . The 
detailed analysis that justifies our assertions is carried out 
in Buchert et al. 



1997 



For equations of state of the form p = up, the time- 
dependent mass term ( |l6| ) is identically zero for flat FRW 
backgrounds, but evolves with time in closed and open 
FRW backgrounds, keeping a constant sign during the 
epoch of interest: m 2 > in open backgrounds and m 2 < 
when the background is closed. 

There are two time scales in our problem: an intrinsic 
time scale associated with the dynamical evolution pre- 
scribed by equation (filf), which is determined by the val- 
ues of the coefficients in it, and an extrinsic time scale 
associated with the background expansion, which is de- 
termined by the parameters defining the background cos- 
mology. A numerical study of the function m 2 (t) for the 
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physically interesting range of parameters, reveals that 
its relative variation over these two time scales is small, 
which justifies an "adiabatic" approximation: we can as- 
sume m 2 (t) to be a constant, rather than a function of 
time, when performing the Renormalization Group analy- 
sis in the next section. This assumption greatly simplifies 
the analysis, and the presence of a non-zero mass term in 
the KPZ equation results in a much richer renormalization 
group trajectory flow and fixed point structure than that 
corresponding to the massless case, as we will see below. 



4. The Renormalization Group Equations 

The starting point for the subsequent analysis is the KPZ 
equation for the velocity potential, eq. (|ll|). Because this 
is a nonlinear equation, and because there are fluctua- 
tions, as represented by the stochastic noise term, it is 
to be expected on general grounds that renormalization 
effects will modify the coefficients appearing in this equa- 
tion. Due to the existence of fluctuations, the coefficients 
in equation ( p"3| ) (y, A, to 2 ) do not remain constant with 
scale: as is well known (Ma |l976| ; Amit |1984[ Binney |1992[ ) 
there are renormalization effects that take place and mod- 
ify these coefficients (or "coupling" constants); the mod- 
ifications are calculable and can be computed using dy- 
namical renormalization group (DynRG) techniques. The 
renormalization group is a standard tool developed for re- 
vealing in a systematic way, how couplings change with 
scale in any physical system under the action of (coarse)- 
graining. Moreover, it predicts that all correlation func- 
tions go as power laws when the system is near a fixed, 
or critical point, i.e., they exhibit self-similar behavior. 
However, the value of the scaling exponent changes from 
one fixed point to another. In fact, the values for the cou- 
plings at some reference scale establish the fixed point to 
which the system will be attracted to or repelled from, be- 
cause these reference values will belong to a specific "basin 
of attraction" . We will look for the IR stable fixed points 
since these reflect the system behavior characteristic of 
long times and large distances and ceases to change as we 
look at the system at ever larger scales. 

Indeed, even before fluctuations are accounted for, one 
can easily obtain the so-called canonical scaling laws for 
the parameters appearing in (|l3|); these follow by per- 
forming the simultaneous scaling transformation x — > sx, 



t —> s z t and ip — > s x i(j and requiring the resultant equa- 
tion to be form-invariant, that is, that the KPZ equation 
transforms to a KPZ equation at the new, larger (s > 1) 
length and time scale. This simple requirement leads to 
scale dependence in the parameters appearing in the orig- 
inal equation of motion as well as in the correlation func- 
tions built up from the KPZ field ip(x,t). When fluctu- 
ations and noise are "turned on" , renormalization effects 
change the scaling behavior of the couplings and correla- 
tion functions away from their canonical forms, in a way 



which depends on the basins of attraction for the fixed 
points of the DynRG equations of the couplings. 

4-1. Calculation of the RGEs. 

The RG transformation consists of an averaging over 
modes with momenta k in the range A/s < k < A where 
s > 1 is the scale factor for the transformation, followed 
by a dilatation of the length scale x — > sx in order to bring 
the system back to its original sizeQ. Here, A plays the role 
of an UV-momentum, or short distance, cut-off character- 
izing the smallest resolvable detail whose physics is to be 
described by the dynamical equations. In classical hydro- 
dynamics, this scale would be identified with the scale in 
the fluid at and below which the molecular granularity 
of the medium becomes manifest and the hydrodynamic 
limit breaks down. The renormalization group transforma- 
tion has the important property that it becomes an exact 
symmetry of the physical system under study whenever 
that system approaches or is near a critical point, because 
it is near the critical point (or points) where the system 
exhibits scale-invariance or, equivalently, where the sys- 
tem can be described in terms of a conformal field theory. 
The hallmark for a system near criticality is that its corre- 
lation functions display power law behavior. This means 
that the statistical properties of the system remain the 
same, except possibly up to a global dilatation or change 
of unit of length. We can calculate the power law expo- 
nents by requiring that the dynamical equations remain 
invariant under the above RG transformation and under 
the further change of scale 



sx 



t — > s z t, and * 



(18) 



where z and \ are numbers which account for the response 
to the re-scaling. By eliminating s from the above, we 
arrive at the fact that the two-point correlation function 
for the velocity potential scales as 

<*(*, *)*(*', 0) «| X - x' | 2 * / ^Zt\z ) > ( 19 ) 

where x is the roughness exponent, z the dynamical ex- 
ponent, and the scaling function f(u) has the following 
asymptotic behavior (see, e. g., Barabasi & Stanley 1995): 



lim f(u) -> u 2x/z , 



lim f(u) — > constant. 

u— *o 



(20) 
(21) 



Notice that because of eq. (fl8|), large s means going 
to the large distance or infrared (IR) and (for z > 0) long 

2 The term "renormalization group" is actually a misnomer, 
since due to the process of averaging over small distances and 
the subsequent loss of information, the RG transformations do 
not have an inverse. Consequently the technically correct name 
would have to be the "renormalization semi-group" . 
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time limit, while small s corresponds to the short distance 
or ultraviolet (UV) and short time limit. 

Making use of the relations Eq.(@) and Eq. @, one 
can immediately obtain the scaling behavior of the 2— 
point correlation function for the density contrast, and 
thus study the asymptotic behavior of this function in 
different regimes. In the following we will implement this 
procedure. First, we obtain the scaling or RG equations 
for the couplings by imposing form-invariance on the equa- 
tion (^), then by using the property of constancy of the 
couplings near fixed points, we obtain and calculate the 
fixed points themselves and the corresponding values of 
the exponents. 

We now characterize the Gaussian noise, r)(x,t). This 
is done by choosing the noise correlation function. Here we 
will use colored or correlated noise, whose Fourier trans- 
form satisfies 



(fj(k,L>)} 

(fj(k^)fj(k',cj')) 
D(k,u) 



0, 

2D(k, Lj)(2ir) d+1 5(k + k')8(uj + u'), 
D Q + D g k- 2p uj- 20 , (22) 



where Dq and Dg are two couplings describing the noise 
strength, and the p, exponents characterize the noise 
power spectrum in the momentum and frequency domains. 
White noise corresponds to Dg = 0. The explicit func- 
tional form of the noise amplitude D(k,iv) reflects the 
fact that correlated noise has power law s ingularities of 
the form written above (Medina el al. 1989 ). 

The solution of ( p"3| ) can be carried out in Fourier space 
where iterative and diagrammatic techniques may be de- 



veloped (Medina et al. 1989). A standard perturbative ex- 
pansion of the solution to Eq. ( |l5| ) coupled with the re- 
quirement of form-invariance and the property of renor- 
malizability, lead to the following RG equations for the 
couplings: 



dm 2 
dlogs 

dv 
dlogs 

dD 
dlogs 

dDg 
dlogs 

dA 
dlogs 



= u[z-2- ^^A d - 2 V- 2 {(d ~ 2V~ 1 )D 
v A Ad 

+(d - 2V- 1 - 2p)D e V~ 2e sec(07r)(l + 29)} 
\ 2 K d 



D (z-2 X -d) + 



A d-2y- 



(23) 



[L>o + 2£> A)F sec(07r)(l + 20) 
+D 2 g V- 40 sec(207r)(l +40)], 

De[z(l + 20)-2 X -d + 2p], 
X[x + z-2 

X 2 Dg Kd^ d -2y-3-2e 



0(1 + 20) sec(7T0)]. 



The full details of the calculation of these equations will 
be presented elsewhere (Martm-Garcfa & Perez-Mercader 
1998| ). The calculation of the RGE's for the massless KPZ 
equation in the presence of colored noise is given in (Med- 
ina et al. 1989). Here, s is the scale factor of Eq. (|l8|), 
and Kd — ^y 1S a dimensionless geometric factor pro- 
portional to the surface area of the d-sphere S d , where d 
is the spatial dimension. This set of equations describes 
how the coupling constants v, A, to 2 , Dq, Dg evolve as one 
varies the scale at which the system is studied, a pro- 
cess commonly referred to in the condensed matter lit- 
erature as "graining" . In the present context we will be 
interested in the coarse-graining behavior, since we seek 
to uncover the behavior of our cosmological fluid as one 
goes to larger and larger scales. In actual practice we will 
investigate the coarse-graining flow in a two-dimensional 
(white noise) or three-dimensional (correlated noise) pa- 
rameter space. This is because one may cast the above 
set of RG equations in terms of a smaller, yet equiva- 
lent set, by employing the dimensionless couplings de- 
fined as V = 1 + ^ , U = \ 2 D K d A d - 2 /iv 3 , and 
U e = \ 2 D e K d A d - 2 - 2 P- i9 /is 3+2 P. Doing so, one arrives at 
the following reduced set of renormalization group equa- 
tions: 



dU 
dlogs 



(2 - d)U 

U Q Ug 



TI 2 
1 o 



[d + 3(dV-2)] 



4dV 3+2e 
[2d- 



4dV 3 
sec(27T0)(l + 40) 

(1 + 20)sec(?r0) x 



' + 3((d - 2p)V - 2)] 



(24) 



dU e 
dlogs 



(2-d + 2p + 40)U g + 



UgUe 
4dV 3 



(dV-2)(3 + 20) 



4dV 3+2f> 



(1 + 20)sec(7r0) x 



dV 
dlogs 



= (v-m+ 



((d-2p)V -2)(3 + 26)} 
1 



: [(dV-2)U 



4dV 3 ' 

-((d - 2p)V - 2)V- 26 {\ + 20) sec(7T0)[/e]] 



The dimensionless couplings Uq , Ug measure the 
strength of the "roughening" effect due to the combined 
action of noise and the non-linearity against the "smooth- 
ing" tendency of the diffusive term in equation ([l3]) : these 
couplings grow when either A or the noise intensity (as 
given by the parameters Dq and Dg) grow, and they be- 
come smaller when the diffusive action (measured by v) 
grows. The dimensionless coupling V measures the com- 
petition between the diffusion and mass terms. Indeed, 
neglecting for the moment the noise and nonlinear terms, 
one can write (fL3|) in Fourier space as 



(25) 



— = -(yk l + TO 2 )*. 
OT 
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Now, if m 2 > (i.e., V > 1), the perturbations in the field 
4" are damped. But if on the other hand, m 2 < (V < 1), 
there exists a length scale L = \J —-j^i below which the 
perturbations are damped, but above which they grow. 
The case V < corresponds to a length scale L < A -1 : 
this scale is smaller than the minimum resolvable length 
scale in the problem at hand, and therefore the diffusive 
term becomes unimportant. We remark in passing that 
the limiting behavior in the RG flow as V — > + leads to a 
technical requirement: one must assume that the diffusive 
term is non-negligible (i.e., v ^ 0) in order to compute the 
RG equations. 

From the above equations one calculates their fixed 
points to which the graining flow drives the couplings. The 
corresponding fixed point exponents z and \ follow from 
substituting the fixed point solutions so obtained back into 
the previous set of RG equations (|||). These exponents 
control and determine the calculated asymptotic behavior 
of the correlation functions. The system will be attracted 
to fixed points in the IR or UV regimes depending on 
whether they are IR-attractive or UV-attractive although, 
as with any autonomous set of differential equations, other 
possibilities exist. 



5. Fixed Points. 

The fixed points, collectively denoted by g* , are the so- 
lutions to the system of algebraic equations obtained by 
setting the right hand sides of equations (pij ) equal to 
zero. They control the stability properties of the system 
of DynRG equations and therefore they characterize the 
asymptotic behavior of correlation functions. When the 
couplings attain their critical values, the system becomes 
critical and the correlation functions enter into power law 
regimes, characteristic of physical systems in the critical 
state where they are scale invariant. Which fixed point the 
system is attracted to or repelled from depends on the basin 
of attraction where the initial conditions for the DynRG 
equations are located. As the couplings evolve under grain- 
ing into their critical values, the n-point correlation func- 
tions approach power law (or scaling) forms because they 
are the solution to first order quasi-linear partial differen- 
tial equations (the Callan-Symanzik (CS) equations (Amit 
1984; Binney [1992 )) whose characteristics are the solu- 
tions to Eqns. (1231). This CS equation simply expresses the 



fact that the correlation functions are invariant under the 
renormalization group transformation. In other words, the 
system is scale-invariant at its fixed points, hence, correla- 
tion functions must take the form of algebraic power laws, 
since these are the only mathematical functions that are 
scale-invariant. The physics of each of these fixed points 
is different, because it depends to a large extent on the 
magnitude and sign of the scaling exponent as well as on 
the attractive or repulsive nature of the fixed point. 

The flows and the fixed points of our set of RG equa- 
tions are represented in Fig. [l], for white noise and Fig. ^| 



V 




Fig. 



1. Flow diagram in (Uq,V) space. Arrows indicate IR 
oo) flow. 




Fig. 2. Flow diagram in (Ug, Uq, V) space around P5+. 



for colored noise, while the features of our fixed point anal- 
ysis are conveniently summarized in the adjoining Table 
1. 

5.1. The value of the Critical Exponents at each Critical 
Point. White Noise. 

In the case of white noise (Dg = 0) we solve for the 
fixed points of the (reduced set) of RG equations. Since 
Ug = in this case, these fixed points will all lie in a two- 
dimensional coupling space spanned by (Uq, V). There are 
three fixed points, labelled as P1,P2 and P3, whose co- 
ordinates are listed in Table 1. Two of them, P1,P2 lie 
in the line V* = 1, corresponding to m* 2 = in the RG- 
evolved KPZ equation, while -P3 lies in the region bounded 
between the V = 1 and V = lines. By substituting these 
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Table 1. Characteristics of the fixed points. 



Point 


Position (Ug,Uo,V*) 


( 2 ,x) 


7 = 4 - 2 X 


IR-eigenvalues 


Class 


PI 


(0,0,1) 


(2,-1/2) 


5 


(4.7,-1,2) 


Saddle point 


P2 


(0,2,1) 


(13/6,-1/6) 


13/3 


(5.2,1,2.2) 


IR-repulsive 


P3 


(0,1.55,0.38) 


(0,2) 





(-1.7,2.0,-15.1) 


Saddle point 


P4± 


V* = 1, [US, US) depend on (p,9) 
Some pairs are not allowed by a dis- 
criminant condition* 


Depend on (p, 9) 
X < or X > 4/3 


Depends on (p, 9) 
7 > 4 or 7 < 4/3 


Sign depends on 






For p = 2.65, 9 = 0.1 












P4+ (4.2,6.4,1) 


(0.64, 1.53) 


0.94 


(-6.2, 1.8,0.64) 


Saddle point 




P4- (3.4, 1.7, 1) 


(0.62,1.52) 


0.96 


(-5.5,-1.6,0.62) 


Saddle point 


P5± 


{Ue,U5,V) depend on (p,9) Some 
pairs are not allowed by a discriminant 
condition*. Essentially the allowed re- 
gion is p e [0,7/2], 
9 G [0, 1/4) 


(0,p-3/2) X G 
[-3/2,2] which in- 
cludes the range 
(1.1, 1.25) for X 


7-2p 


Sign depends on 






For p = 2.65, 9 = 0.1 












P5+ (-0.1,0.43,0.19) 


(0,1.15) 


1.7 


(-0.4U3.5i, -27.3) 


IR-attractive 




P5- (3.0,-1.96,-0.55) 


(0,1.15) 


1.7 


(-19.2,8.5,-1.5) 


Saddle point 



* These fixed points are obtained as real solutions of a second degree algebraic equation, and therefore the set of pairs that give 

real fixed points are restricted by a discriminant inequality. 

Remarks: 



— Only P5+ has a basin of attraction, which is approximately in the region V £ (0, 1), Uo € (0, 1.5), Ug £ (—0.2, 0). P2 has a 
large "basin of repulsion", and the other points are saddle points. 

— The trajectories outside the P5+ basin of attraction reach infinity or fall to the singular V = plane at a finite "time" . 

— Only with P5+ can we adjust the observed values 7 £ (1.5, 1.8), because P5- has V* < in this region, which implies 
smaller length scales than the one corresponding to the UV cut-off A. 



fixed point values back into the original set of RG equa- 



tions (23), we solve for the corresponding fixed point expo- 
nents x an d z. These values are listed in the third column 
of Table 1 . Linearizing the RG equations about each one 
of their fixed points allows one to calculate the infrared 
(IR) stability properties of the fixed points and thus char- 
acterize their behavior with respect to the coarse-graining. 
This entails expanding the RG equations in C/ — > Uq+SUq 
and Vq — > Vq + SVq to first order in the perturbations and 
solving for 5U and SV . In general, the eigenvectors of 
this linearized system will involve linear combinations of 
SUq and SVq. The associated eigenvalues, which determine 
the stability properties of the fixed point, are also listed in 
Table 1. A positive eigenvalue indicates that the coarse- 
graining induces a flow away from the point along the 
eigen-direction, while a negative eigenvalue indicates the 
point is stable in the infrared, since the flow will be into the 
fixed point. The type, or class, of fixed point, whether it be 
IR-attractive (all eigenvalues negative), IR-repulsive (all 
eigenvalues positive) or a saddle point (mixed sign eigen- 
values) is listed in Table 1. The coarse-graining flow of the 
couplings in the white noise case is given in Fig. . These 
flow lines were calculated by integrating numerically the 
set of two coupled first order differential equations for Uq 
and V, using general choices for the RG equation initial 
conditions. 



5.2. The value of the Critical Exponents at each Critical 
Point. Colored Noise. 

For colored or correlated noise (De ^ 0), we now deal 
with a three dimensional space of dimensionless couplings 
Ug,Uo and V. Using the same procedure as discussed 
above, we solve for the RG equation fixed points which 
lead to a total of seven fixed points, including the same 
three points P1,P2 and P3 obtained in the white noise 
limit. Thus, allowing for colored noise yields four addi- 
tional fixed points which we will denote by P4± and P5±, 
since they arise in pairs. Their coordinates are listed in the 
Table, together with their associated exponents x an d z 
and IR eigenvalues and stability properties under coarse- 
graining. Both, the positions of these points and values 
of their exponents, depend in general upon the values of 
the noise exponents p and 9, which parametrize spatial 
and temporal correlations in the noise. The pair of fixed 
points labelled as P4± always lie in the plane V* = 1, 
(i.e., m* 2 — 0) but their location within this plane varies 
with the noise exponents. Moreover, any RG flow which 
starts off in this plane will always remain in this plane 
(this plane acts as a separatrix). It is important to point 
out that there are values of the noise exponents p, 9 which 
lead to complex values of the fixed point coordinates for 
P4±. This is because these fixed points arise as solutions 
of a quadratic algebraic equation whose discriminant can 
become negative for values of p and in certain domains 
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in parameter space. We must rule out such values of the 
noise exponents because of physical reasons. For the al- 
lowed values of p and 8 (i.e., those that lead to real fixed 
points) , wc find that P4± have exponents z and x that are 
complicated functions of these parameters, but we have 
checked that for all allowed values, either \ < or \ > f ■ 
The IR eigenvalues and the nature of these fixed points 
depend on p, 8. For the choice shown, i.e., p = 2.65 and 
8 = 0.1, both P4± are saddle points. In fact, we have 
confirmed that the pair P4± will always be saddle points 
whenever p € (2.60, 2.75) and 8 € (0, .23). The reason for 
choosing these particular exponent intervals will become 
clear in Section || . 

The next pair of fixed points P5± also have coordi- 
nates and exponents depending on the noise exponents. 
As in the case of P4±, there are values of p, 8 (same as 
for P4±) leading to complex couplings, which we rule out. 
The allowed values of the noise exponents lie roughly in 
the respective intervals p G [0, |] and 8 € [0, j). The ex- 
ponents for P5± are given by the simple functions z = 
and x = P ~ | ! an d the stability properties do not depend 
on the particular values of the noise exponents within the 
above mentioned intervals. For the choice p = 2.65 and 
8 = 0.1, P5+ is IR attractive while P5— is an unstable 
saddle point. We must nonetheless exclude P5— from our 
consideration since it lies below the plane V — 0, i.e., it 
corresponds to a renormalized KPZ mass m* 2 < —vA 2 , 
which therefore exceeds the scale of the momentum cut- 
off imposed on our perturbative calculation. 

6. Scaling properties of matter density 
correlations. 

We can now begin to put together the results we have 
obtained in relation with the scaling properties predicted 
by the DynRG. Using equations (|l^) and (lj) we can write 
for the 2-point correlation function of the density contrast 
in comoving coordinates that 



which on using equation (Eq) above for equal times gives 



£(r) 



for the correlation function, and 



P{k) ~ k 2+ \ 



(30) 



(31) 



for its Fourier transform, the power spectrum. 

Equations ( j30|) and_(31) are now suitable for compar- 



ison with observations^, since they are written in comov- 
ing coordinates, and therefore correspond to quantities di- 
rectly inferred from observations and provided we assume 
that light traces mass. Otherwise we would need to cor- 
rect for bias, a strategy we leave for a future publication. 
From d30| ) and the definition of the exponent rj, we see that 
the exponent 7 measured from large scale galaxy surveys, 
where £,OBs( r ) ~ r ~ 7 , is calculable in term of the rough- 
ening exponent x an( i is given by Q 



7 = 4-2*(p,0). 



(32) 



Using our fixed point analysis we have computed all 
the exponents x an d z for all the fixed points (fourth col- 
umn in Table 1). The predicted values for 7 are tabu- 
lated in the Table as shown. Thus, at each fixed point, 
we derive a power law for the density correlation function 
f(r) ~ r~ 7 . In the case of white noise alone, we see that 
none of the three corresponding fixed points are capable of 
reproducing any of the inferred values of 7 (1.5 ~ 7 ~ 1.8) 
from observations. It is of interest to point out, however, 
that the point P3 yields the exponent 7 = 0, which means 
that the correlation function is strictly constant at this 
point, i.e., the matter distribution is perfectly homoge- 
neous. However, P3 is an IR unstable saddle point and 
a fine-tuning in the initial values of Uo and V would be 
required in order to have the system flow into it under 
coarse- graining. 



a\x-y\,t) = (S(x,t)S(y,t)) 
( F \ 2 f 2 



(26) 



We see that the scaling behavior of £(r) is determined once 
we know the scaling behavior for the equal time two-point 
correlation function of the density contrast for the velocity 
potential ^(x, t), 

C(\x-y\) = (*(x,tMv,t)) . (27) 

iFrom Eqns. @ and @, we have 

(9(x,t)9(y,t)) ~ \x-y\ 2 *. (28) 

Introducing (for convenience and notational consistency 
with the literature in critical phenomena) 77 = 2 — d — 1\, 
we have 



C(r) 



„-(d-2-H7) 



C(k) ~ k r >- 



(29) 



3 Observation has shown that large scale structure depends 
on cosmological parameters such as the age of the Universe to 
and initial density f2o- While these parameters (and others) 
'do not appear explicitly in our calculations, they are there 
implicitely in quantities such as ro, the scale factor, whose 
calculation cannot be carried out without a deeper knowledge 
of the noise sources. Nevertheless, the renormalization group 
gives one a way to calculate the power law scaling law and its 
exponent as a function of which particular basin of attraction 
the systems happens to start out in. This explains exactly what 
we mean when we claim that no fine tuning of parameters 
is required in order to account for the value for the scaling 
exponent inferred from the observations. Irrespective of where 
the system starts out in a particular basin of attraction, in 
time it inevitably flows into that basin's fixed point: this is 
very much a kind of independence of initial conditions. 

4 Notice that when written in this way, the value of 7 seems 
to be independent of the spatial dimensionality d. But this is 
fallacious, since x itself depends on d as is seen from the RGE's. 
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For colored noise, the RG behavior becomes signifi- 
cantly richer since we now must map out RG equation 
flow trajectories in a three dimensional space of couplings 
(see Fig. ||). A glance at the tabulated calculated val- 
ues of the exponent 7 shows that P5+ is the only fixed 
point capable of yielding values of 7 within the currently 
accepted range of values inferred from observations, by 
choosing p E (2.60,2.75) and any 9 E (0, .23), since 
X — x(p) = 7 — 2p for P5±. Although 7 is independent of 

6, a non-zero value of 9 in the above interval must be cho- 
sen in order that P5± exist (this pair of fixed points vanish 
identically when 9 = 0) . It is also the only IR-stable fixed 
point in the three-dimensional coupling space spanned by 
the couplings (V, Uq, Ug). This implies the following be- 
havior: any point initially in the basin of attraction of P5+ 
inevitably ends up at this fixed point as larger and larger 
scales are probed: no fine tuning is required. 

7. Discussion and Conclusions. 

In this paper we have presented an analytical calcula- 
tion of the density-density correlation function for non- 
relativistic matter in a FRW background cosmology. The 
calculation hinges on blending two essential theoretical 
frameworks as input. First, we take and then reduce the 
complete set of non-relativistic hydrodynamic equations 
(Euler, Continuity and Poisson) for a self-gravitating fluid 
to a single stochastic equation for the peculiar- velocity po- 
tential tp. The noise term rj(x,t) models in a phenomeno- 
logical but powerful way different stochastic processes on 
various length and time scales. Second, we apply the well- 
established techniques of the dynamical Renormalization 
Group to calculate the long-time, long-distance behavior 
of the correlation function of the velocity potential (which 
relates directly, as we have seen, and under the conditions 
we have specified, to the density-density correlation) as a 
function of the stochastic noise source in the cosmological 
KPZ equation. 

When we carry out a detailed RG fixed point analy- 
sis for our KPZ equation, we find that simple white noise 
alone is not sufficient to account for the scaling exponent 
(1.6 ~ 7 ~ 1.8) inferred from observations of the galaxy- 
galaxy correlation function. However, we can get close to 
this range for the observed exponent if the cosmic hydro- 
dynamics is driven by correlated noise. In fact, we need 
only adjust the degree of spatial correlation of the noise 
to achieve this, since the calculated value of the exponent 
7 = 4 — 2x(p) associated with P5+ is independent of the 
degree of temporal correlations, as encoded in 9. More- 
over, P5+ is the only fixed point with this property, and 
is simultaneously IR-attractive, both desirable properties 
from the phenomenological point of view. This last prop- 
erty is very important, since it means that the self-similar 
behavior of correlations is a generic outcome of the dy- 
namical evolution, rather than an atypical property that 
the system exhibits only under very special conditions. 



Thus we have learned that colored noise seems to play 
an important role in the statistics of large scale structure. 
Now, why it is that a non-zero 9 and a p E (2.6, 2.75) 
are the relevant noise exponents is a question one still 
has to ask. To answer it, one would in principle have to 
derive the noise source itself starting from the relevant 
physics responsible for generating the fluctuations. Here, 
we content ourselves with the phenomenological approach. 
Nevertheless, at some future point, it may be possible to 
say more about the spectrum of noise fluctuations. 

It must be pointed out that the Renormalization 
Group allows us to compute only asymptotic correlations, 
i.e., after sufficient time has elapsed so that the effect of 
noise completely washes out any traces of the initial condi- 
tions in the velocity field or initial density perturbations. 
Hence, there will be a time-dependent maximum length 
scale, above which noise has not become dominant yet and 
correlations therefore remember the initial conditions, as 
well as a transition region between these two regimes. The 
length scale at which this transition occurs depends on 
noise intensity and on the amplitude of initial perturba- 
tions. Therefore, a physical model of the origin of noise 
would enable us to predict this scale, but such a task is 
well beyond the scope of the present paper. 

It is also interesting that the Renormalization Group 
allows us to compute the proportionality coefficient of 
the matter density correlations, i.e., the quantity ro in 
£(r) = (r/Vo) -7 , by means of the so-called improved per- 
turbation expansion (Weinberg 1996| ). This length scale 
ro marks the transition from inhomogeneity (r <c ro, 
£(r) 3> 1) to homogeneity (r 3> ro, £(r) <C 1). But r 
is also dependent on the noise intensity and other free 
parameters of our model, so that a prediction is impossi- 
ble without a deeper knowledge of the noise sources. At 
this point, we should mention the recent debate about the 
observational value of ro: on the one hand there is the 
viewpoint that homogeneity has been reached well within 
the maximum scale reached by observations (Davis 1996) 
(ro ~ bOMpc); on the other hand there is an opposing 
viewpoint as expressed by Labini et al. 1998 who argue 



that ro ~ 200Mpc. Either viewpoint can be easily incor- 
porated into our model, by suitably choosing the noise in- 
tensity. But if the latter value of ro turns out to be correct, 
this would indicate within our model that the noise inten- 
sity is really large and therefore, that noise and stochastic 
processes played a much more important role in the His- 
tory of the Universe than has been thought to date. 
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